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Abstract 

We investigate the behaviour of the lowest nonhydrodynamic modes in a class 
of holographic models which exhibit an equation of state closely mimicking the 
one determined from lattice QCD. We compute the lowest quasinormal mode 
frequencies for a range of scalar self-interaction potentials and find that the 
damping of the quasinormal modes at the phase transition/crossover falls off by 
a factor of around two from conformality after factoring out standard conformal 
temperature dependence. The damping encoded in the imaginary part of the 
frequencies turns out to be correlated with the speed of sound and is basically 
independent of the UV details of the model. We also find that the dynamics 
of the nonhydrodynamic degrees of freedom remains ultralocal, even to a higher 
degree, as we deviate from conformality. These results indicate that the role of 
nonhydrodynamic degrees of freedom in the vicinity of the crossover transition 
may be enhanced. 
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1 Introduction 

The quark-gluon plasma produced in relativistic heavy-ion collisions at RHIC and 
LHC is very successfully described by phenomenological hydrodynamic models [1]. 
Nevertheless, it is quite clear that the hydrodynamic description is applicable only from 
a certain initial time after an inherently nonequilibrium initial phase of expansion. The 
estimation of this so-called ‘thermalization time’ has been at the focal point of much 
theoretical effort. 

The main fundamental difficulty in studying such a problem from first principles is 
that when the plasma is strongly coupled we need a method which would be at the 
same time nonperturbative and which would work directly in Minkowski signature. 
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For these reasons, the methods of the AdS/CFT correspondence have been applied in 
this context although predominantly for the case of the conformal A7 = 4 SYM theory 
[2- 3j. Despite that, the results obtained within this framework are very encouraging. 

Firstly, the appearance of a hydro dynamical description of the resulting plasma system 
is not an input but rather a result from a much more general dual gravitational de¬ 
scription which incorporates genuine nonhydrodynamical degrees of freedom. Thus one 
can study the transition to hydrodynamics and its properties. Secondly, it has turned 
out that at the transition to hydrodynamics, the components of the energy-momentum 
tensor in the local rest frame are still significantly anisotropic, meaning that the plasma 
is then still quite far from thermal equilibrium which indicates that the phrase ‘early 
thermalization’ used in this context is really a misnomer. Thirdly, for numerous initial 
conditions, the plasma behaves hydrodynamically with very good accuracy when the 
dimensionless product of the proper time and temperature Tt ~ 0.6 — 0. 7M- Before 
that nonequilibrium degrees of freedom^ are typically very important. 

It is worth mentioning also numerous important results on the dynamics of shock wave 
collisions pen ei which we do not describe in more detail here. 

Since the full nonlinear dynamics in the deeply nonequilibrium regime is very compli¬ 
cated as it is described by higher dimensional Einstein’s equations and can be studied 
essentially only using the methods of numerical general relativity, it was suggested in 
0HQ] that it may be useful to incorporate just the lowest, least damped nonhydrody¬ 
namic degrees of freedom into the commonly used nonlinear hydrodynamic description. 
On the dual gravity side these degrees of freedom are the so-called quasinormal modes 
(QNM) of a finite temperature black hole. The 4D equations involving these nonequi¬ 
librium modes proposed in mi take as an input from the gravitational description 
only their real and imaginary frequencies coqnm = + iu>i- Moreover, it turns out 

that the dependence of these frequencies on the momentum k is very mild and can be 
neglected in a first approximation. This property leads to a certain ‘ultralocality’ of 
the dynamics of the nonequilibrium modes on top of a hydrodynamic flow. 

All the above investigations were performed in the context of the conformal J\f = 4 
SYM theory which, by its very definition does not exhibit any kind of phase transi¬ 
tion or crossover behaviour. It is thus very interesting to study what modifications 
to the above picture appear for a nonconformal theory. There are various ways to 

1 In the present paper we use this term to denote all nonhydrodynamical degrees of freedom in the 
plasma. 
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model nonconformal theories within the AdS/CFT correspondence either following a 
top-down approach by studying a specific nonconformal theory with a known string 
theory construction, or a bottom-up approach where the gravitational background is 
phenomenologically fixed to give properties known from lattice QCD. In the present 
paper we decided to concentrate on the latter approach partly for simplicity and partly 
in order to deal with a gravitational system which has very similar equation of state 
to real QCD^j 

In this paper we will concentrate on the dynamics of the lowest nonhydrodynamic 
degrees of freedom in the nonconformal setting. In particular we will investigate how 
the damping of these modes changes when we approach the crossover temperature T c 
(defined more precisely later). This answers an important question whether the role 
of nonhydro dynamic degrees of freedom becomes more important or less important 
closer to the phase transition. Secondly, we will investigate whether the ‘ultralocality’ 
property observed for AT = 4 SYM nonequilibrium degrees of freedom (QNM) still 
holds in the nonconformal case, especially close to the crossover/phase transition. 

The plan of this paper is as follows. In section 2 we will review the family of gravita¬ 
tional backgrounds that we will consider. In sections 3 and 4 we will give some details 
on their explicit numerical construction and on our procedure for finding quasinormal 
mode frequencies. As far as we know such a procedure has not been employed so far 
in the literature and may be useful also in other contexts. There we will also discuss 
the relation of scalar QNM with the metric ones for these backgrounds and proceed, 
in section 5, to describe our results. We close the paper with conclusions. 


2 The nonconformal gravity backgrounds 


In this paper we will study a family of black hole backgrounds which follow from an 
action of gravity coupled to a single scalar held with a specific self interaction potential: 


S 



R~\ W ) 2 + V(<P) 


( 1 ) 


This family has been introduced by Gubser and collaborators in a series of papers 
mmm and used to mimic the QCD equation of state by a judicious choice of the 
potential. These backgrounds were then used to study bulk viscosity (which identically 

2 A paper investigating the complementary top-down approach appeared simultaneously 12 a. 
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Figure 1: The speed of sound c 2 for the potentials Vi — V 5 correspond to deformations 
of the theory by operators of dimensions 3.93, 3.67, 3.55, 3.10, 3.00 respectively given 
in table 1, together with Lattice QCD data from j!6j. 


vanishes in the conformal case) and quite recently in (I3J, where second order hydrody¬ 
namic transport coefficients have been calculated. Ref. D2S provided an approximate 
but quite accurate formula for computing the equation of state (or more precisely the 
speed of sound c 2 ), 


c 2 _ d log T 


d log S' ^ 

directly in terms of the scalar potential V(</>). Thus the scalar potential parametrizes 
the physics of the particular type of gauge theory plasma. 


The family of scalar potentials that we consider is 


V (0) = —12 cosh (6 0) + C 2 0 2 + C 4 0 4 + Cg 0 6 (3) 


The quadratic terms in 0 (mass term) determine, according to the standard AdS/CFT 
dictionary, the dimension A of the operator A (A — 4) = m 2 L 2 , where L is the AdS 
radius which we fix it to one. 

I 11 this paper we will mostly concentrate on a set of parameter choices in ([3]) which 
approximately reproduce the equation of state of QCD plasma as determined by the 
Budapest-Wuppertal group in [16]. The resulting potentials are listed in table [l] and 
the corresponding speed of sound is shown in Fig. [l] In that figure we show the 
speed of sound extracted from a numerical construction of the corresponding black 
hole solution, as described in the following section, together with the lattice QCD data 
for c 2 . Note that in each case we are free to choose the units of temperature. Here, 
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potential 

b 

c 2 

c 4 

c 6 

A 

Vx 

0.606 

2.06 

-0.1 

0.0034 

3.93 

V2 

0.606 

1.6 

-0.1 

0.0034 

3.67 

Vz 

0.606 

1.4 

-0.1 

0.0034 

3.55 

v 4 

0.606 

0.808 

-0.1 

0.0034 

3.10 

V 5 

0.606 

0.703 

-0.1 

0.0034 

3.00 


Table 1: The first group of potentials mimicking the QCD equation of state, together 
with the dimensions A of the relevant scalar operators. Potential V fi in table 2 was 
first constructed in [T2], while V 5 was recently used in [15]. 


following ra we fix this freedom so that the temperature corresponding to the lowest 
dip in c 2 coincides for all the potentials. This will also be our provisional definition 
of the critical temperature T c . According to [15] this value should be 143.8 MeV for 
QCD. Note finally that for high temperatures, the equation of state becomes essentially 
conformal. Similarly other properties such as the quasinormal frequencies also approach 
the conformal values characteristic of AT = 4 SYM at high temperature. 

In order to check that the qualitative conclusions are generic, we also considered some 
other potentials leading to different profiles of c 2 (T). We will discuss them in sec¬ 
tion 15.41 


3 Eddington-Finkelstein coordinates 

This section describes the black hole background solutions for the quasinormal mode 
calculations. These backgrounds are the same as those in [12], but since our goal is 
to determine the quasinormal mode frequencies, it will be convenient to express them 
in Eddington-Finkelstein coordinates, rather than in the coordinates used in |T2|. We 
will discuss this in more detail in the following section on quasinormal modes. 

The Ansatz for these solutions follows from the assumed symmetries: translation in¬ 
variance in the Minkowski directions as well as SO{3) rotation symmetry in the spatial 
part. This leads to the following form of the line element: 

ds 2 = g u dt 2 + g xx dx 2 + g rr dr 2 + 2g rt drdt (4) 

where all the metric coefficients appearing in Q are functions of the radial coordinate 
r alone, as is the scalar field 0. This form of the field Ansatz (determined so far only by 
the assumed symmetries) allows two gauge choices to be made. In JI2] a Schwarzschild- 
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like gauge was adopted by taking g tr = 0, accompanied by the condition 0 = r. The 
latter condition on the scalar field lead to key simplifications which were used to solve 
the field equations. The final form of the Ansatz in H3 was thus 


ds 2 = e 2A (—hdt 2 + Ax 2 ) H——dr 2 (5) 

h 

(f) = r (6) 


where A, B , and h are functions of r (or, equivalently 0). 

For the purpose of computing the quasinormal modes it is very convenient to use a 
different gauge - the Eddington-Finkelstein gauge g rr = 0. It is typically convenient 
to also impose the gauge choice g tr = 1 , but for our purposes it turns out to be very 
effective to use the remaining gauge freedom to set 0 = r. Furthermore, if we label the 
metric components as 


ds 2 = e 2A (—hdt 2 + Ax 2 ) — 2 e A+B dtdr (7) 

0 = r (8) 


then the field equations take the form 

A" - A'B' + - = 0 (9) 

6 

h" + (4 A 1 - B')ti = 0 (10) 

6 A'ti + h(24A' 2 - 1) + 2e 2B V = 0 (11) 

h' e 2B 

4 A'-B' + -- —V' = 0, (12) 

h h 


where the prime denotes a derivative with respect to 0. 

With the assumed labelling of metric coefficients ([5]), equations ([£])-( 12) are identical to 
those of appearing in [12]. and so they can be solved following the method described 
there (see also USD- In the remainder of this section we review this procedure for 
completeness. 


We are interested in solutions possessing a horizon, which requires that the function h 
should have a zero at some 0 = 0#: 


K^h) = 0, (13) 

It is easy to see that the solutions of Eq. 0 -([l2| can be expressed in terms of a single 
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function G(0) = A'(eft): 


= B 


H 


A{4>) =A h + d 0G(0), 
J (pH 

G(0) \ t [* d0 


In 


+ 


G((Ph)J )<p H 6G(0) ’ 


/i(0) = h H + hi / d0e- 4A ^ )+s( ^ 

J (pH 


(14) 

(15) 

(16) 


In the expressions above Ah, Bh, hn and h\ denote constants of integration which will 


be determined by requiring the appropriate near-boundary behaviour and eq. (13). 


As in [12] by manipulating the held equations (|9])- (12) one finds the nonlinear “master 
equation 1 ' 

rj H fry i rj \ 

(17) 


G’ _ d (G 1 G’ 

G + V/W ~ d0 ln V G + 6G G + V/W 

The strategy is to solve this equation numerically by integrating it from the horizon at 
0 = (j)H down toward the boundary at 0 = 0. Once G is known, the metric coefficients 


can be recovered from Eq. (14)-(16). 


Solving Eq. ( flT| ) requires appropriate boundary conditions, which can be determined 
by evaluating ([IT]) and (12) at the horizon and using (13). In this way one finds 


V(cPh) = —3 e- 2S ^G(0^)/0(0i0, E'(0) = e~ 2B ^h'(cp H ). 


From this it follows that 


g(4h) = 


V (0 


H 


315(0 


H 


(18) 


(19) 


Using (19) and (17) one finds the following near-horizon expansion: 


gw reyy ->i m 


3 V'(<j> H ) 6 V V'(<j> H ) 


In particular, the expansion (|20|) implies that 

V(0^)U // (0h) 


G'(M = l (’ 


V'(M 


( 21 ) 


To summarize: to find a numerical solution of 0 we can specify a value for <p H 
and then use the conditions (19) and (21) as boundary conditions for integrating 0- 
There is however one technical complication in performing the numerical integration 





















outlined above: Eq. (19) implies that at the horizon 


G(0„) + V(<l> H )/(3V'(<t, H )) = 0 


( 22 ) 


which makes some terms of (17) singular. Even though such superficially singular 


terms cancel, their presence makes numerical computations troublesome. In order to 
circumvent this difficulty, instead of (pH one can initialize the integration at a point just 


outside the horizon, at 0 = <pn — £h, where en <C 1. Then using (21) one can calculate 


G((pH ~ €h) and G\(pH — €h) and then use these values as the boundary conditions. 
One also needs to regularize at the boundary (0 = 0) by integrating down to a small, 
but finite value <p = e#. 


Having determined G , one can End the metric from (14)—(16). The constants of inte¬ 
gration can be determined following [12]. The result is 

r4>H 



In (pH 

H 

A -4 

•H 

= In (— 

H 

= 0 


Ah = 

B 
h 

hi = 


+ 


d <p 


G((p) — 


4 V{cp 


h 


E(0 )V'{<p H )L 

1 


d(pe~ 4A ^ +B ^ 


+ 


(A - 4)0 
d 0 


60(0) 


(23) 


This way the metric is determined for any given choice of 0#. 

The Beckenstein-Hawking formula for entropy leads to the following expression for the 
entropy density 

s = (24) 


ni 


and the standard argument requiring non-singularity of the Euclidean continuation at 
the horizon gives 

“ A H - B H \ti{r E )\ 


T = - - 

These equations lead to the formula 

2 d log T / dcpH 

C s 


47T 

1 

_ 

~ 3 


1 V\cp H f 


(25) 


(26) 


d\ogs/d(pH 3 2 V(<Ph) 2 

where the latter equality, proposed in [12], is only approximate but works surprisingly 
well. In the present paper when determining the speed of sound c^, we always use 
the exact formula and determine it from d log T/d log s with temperature and entropy 
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extracted from the exact numerical solutions. 


4 Quasinormal modes 

4.1 Introductory remarks 

It is convenient and enlightening to formulate the problem of finding quasinormal 
modes in terms of gauge invariant variables, which are diffeomorphism invariant linear 
combinations of the perturbations. This approach is well known in general relativity, 
and has been adopted in the holographic context in ra, where the conformal case of 
A7 = 4 supersymmetric Yang-Mills theory was considered. The generalisation to the 
non-conformal cases was undertaken in [TO] . In this section we briefly summarise our 
findings in the context of the models under consideration in this paper. 

Under an infinitesimal diffeomorphism transformation, the metric and the scalar held 
fluctuations transform as the metric and the scalar held itself, i.e. 

9iiu —>- 9nv ~~ ^ — V v (27) 

By examining linear combinations of the linearized perturbations one hnds hve gauge 
invariant channels: two shear channels, one scalar channel, one sound channel and one 
bulk channe|^] [T9]. 

We assume the plane wave e ~ lU}t + lkx dependence of the huctuations on boundary coor¬ 
dinates t, x and some nontrivial dependence on the radial coordinate r. In general, the 
equations for the shear and scalar channels are decoupled from the rest, which have the 
same form as the conformal case. However the sound and bulk channels have coupled 
second order equations [T9] . 

In the zero momentum limit, k —> 0, the equation for the sound mode becomes decou¬ 
pled from the bulk mode and at the same time the equations for the other channels 
reduce to the equation for the QNM’s of a massless scalar held, except for the equation 
of the bulk channel which is still coupled to the sound mode. Interestingly, this has 
two advantages. At k = 0 it is enough to find the QNM’s for an external massless 

3 The bulk channel is a linear combination of transverse metric fluctuations and massive scalar 
fluctuation. One reason we call it ’’bulk” channel is it leads to non-zero bulk viscosity as well [Si- 
One could refer to it as the ”non-conformal” channel and to the rest, which are already known in 
conformal case, as conformal channels. 
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scalar field for the conformal channels. On top of that, the coupling of the bulk mode 
with the sound mode means that the former has the same frequency as the latter. So 
in this limit, all the information about the QNM’s of the theory is summarized in the 
QNM’s of an external massless scalar held. 

In view of this, the following analysis is focused on the QNM’s of an external massless 
scalar held in the nonconformal backgrounds under consideration. Of course, in general 
it would still be interesting to study the metric and massive scalar held perturbations 
in detail. 


4.2 QNM of a massless scalar field 


In view of the results described in the previous section, we turn to exploring the effects 
of conformal symmetry breaking by considering the QNM of a massless scalar held T in 
the background ([7]). As discussed earlier, the equation obtained for this case contains 
all the essential elements for QNM perturbations of the background. 

The held equation for a massless scalar is simply the wave equation 

V. 4 V a T = 0 (28) 

Quasinormal modes are solutions of the form 

T = e~ iuit+ikx ^((l)) (29) 


which satisfy the ingoing condition at the horizon, which in the Eddington-Finkelstein 


coordinate system reduces to regularity there. Substituting (29) into Eq. (28) and using 
the form of the background metric given in Eq. ([Tj) leads to the following equation for 
the amplitude ip ((ft): 


(3 G V' + V)tp" - 3e~ A - B G' ( e A+B V' + 2ioj) o' + 3e~ 2A - B G' ( k 2 e B - 3 iuie A G) ip = 0. 

(30) 

where primes denote derivatives with respect to (ft. Note that only the functions A 
and B appear here - the function h drops out. Since we are imposing two boundary 
conditions, this is an eigenvalue problem which can be solved only for specihc values 
of the complex QNM frequency u. 
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4.3 Numerical approach 


Quite generally, the chief advantage of using the Eddington-Finkelstein coordinate sys¬ 
tem for finding quasinormal mode frequencies is twofold. Firstly, the ingoing boundary 
condition at the horizon gets translated just to ordinary regularity of the solution at 
the horizon. Secondly, due to the special form of the temporal part of Eddington- 
Finkelstein metric, the dependence on the mode frequency of the relevant differential 
equation is linear. Hence the problem of finding the quasinormal frequencies amounts 
to solving a linear ODE of the form 


L iT = caLoT 


(31) 


where L\ and L 2 are specific differential operators, with T satisfying essentially Dirich- 
let boundary conditions at the boundary and being regular at the horizon. While 
many approaches to the problem of finding quasinormal modes have been described in 
the literature ca, we believe that the approach we describe here is very effective in 
conjunction with the spectral representation in terms of Chebyshev polynomials 


This representation reduces the task of solving Eq. (31) to a set of linear equations. 


The differential operators appearing in (31) are represented as matrices, and due to the 


linear dependence of this equation on u this reduces to a generalized matrix eigenvalue 
problem which can be solved very efficiently. 


In the case we are studying, the relevant equation is Eq. (30). This equation is indeed 


linear in c 0 . We have implemented the strategy outlined in the previous paragraph 
and verified the stability of the resulting solutions when varying the number of grid 
points in the Chebyshev discretization. The results of these numerical calculations are 
presented in section [5j 


5 Results 


We are interested in the dependence of the QNM frequencies on temperature and k 

we 


in the vicinity of k = 0. Therefore, in line with the discussion in section 4.1 


focus on the quasinormal modes of an external massless scalar held. By use of the 
term “external” we wish to emphasize that this scalar is distinct from the scalar </> 
appearing in the background geometry, as QNM modes of the latter are mixed with 
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Figure 2: The imaginary parts of the lowest quasinormal mode at k = 0 for the 
potentials from table 1 (left). The imaginary part for potential Vi together with the 
“phenomenological” according to Eq. (33) (right). 


the metric perturbation^} 


5.1 The imaginary part of the QNM frequencies — damping 


In the left panel of figure [2] we show the imaginary parts of the QNM frequencies in 
units of temperature which is the natural scale in the problem i.e. 


Irn u! 
2?rT 


(32) 


We observe that the damping significantly decreases (by a factor of 2) close to the 
transition. This shows that in the nonconformal case nonequilibrium dynamics become 
more important close to T c . Moreover we find that the plots basically lie on top of each 
other for the various potentials from table 1. This indicates that the QNM frequencies 
are not sensitive to the fine details of the potentials but are essentially dependent just 
on the equation of state (speed of sound c^(T)), which was the common denominator 
of all the potentials from table 1 . 


In order to parameterize the dependence of the damping on deviation from conformality, 
we propose a phenomenological formula expressing this as a linear combination of (? s — | 
and T-j^c 2 s {T). Specifically, we posit 


Im ui — Im U7 COI1 f 
2 nT 



j)+7'4ar) 


(33) 


where 7 , 7 ' are phenomenological parameters and 


4 See the discussion in section 


4.1 


— 1.373 is the conformal 
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Figure 3: The real part of the lowest quasinormal mode at k = 0 for the potentials 
from table 1 (left) and the ratio of the real to the imaginary part (right). 


limit value. These parameters can be fitted to the numerically calculated difference of 
the damping w.r.t the conformal case. For the potential V 2 in table 1 we got 


7 = -3.729 , 7 ' = 0.452 


(34) 


In figure fright) we show a plot of (linen — Irn ay on f) / (27 tT) together with the fit. This 
two parameter fit is surprisingly good and may be thus used phenomenologically to 
estimate the damping in a nonconformal theory with the QCD equation of state. Since 
the quasinormal frequencies for the family of potentials we used in the left panel of 
figure [2] basically coincide, the single choice of parameters given in Eq. (34) works well 
for all of them. 


5.2 The real part of the QNM frequencies 

It is worth noting that similarly as for the imaginary part of the quasinormal fre¬ 
quencies, the real parts of the frequencies corresponding to the various potentials from 
table 1 are also very close to each other (see figure Ii eft)). This signihes that the 
QNM frequencies are basically insensitive to differences in the UV (since the various 
potentials correspond to different A’s) and are governed by IR physics i.e. essentially 
the equation of state. 

In figure fright), we also see an enhancement of the real part of the frequencies slightly 
below T c . 
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Figure 4: The curvature at k = 0 of the damping frequency, overlaid with a ~ 1.114 
and f3 — 0.342 for the potential Vi in table 1. 

5.3 Ultralocality 

As indicated in the introduction, an interesting property of the dispersion relation for 
the nonhydrodynamic degrees of freedom in the conformal case of A/" = 4 SYM theory 
is the very mild dependence of the frequencies ojj and ojr on the momentum k. In this 
section we show that for nonconformal theories this property holds to an even higher 
degree. Interestingly, the curvature of the damping i.e. 

Imc jj"(k = 0) 


2v tT 


(35) 


follows to a surprising accuracy (at least until T ~ T c , then it starts to deviate) just 
the speed of sound c^, up to an overall numerical factor determined in the conformal 
high temperature limit (i.e. equivalently in J\f = 4 SYM). The relevant plot is shown 
in figure [4} The plot for the other potentials in the same family are basically the same. 


5.4 Further examples 

As discussed in section[2j the choice of scalar potential translates into a specific equation 
of state for the QCD-like gauge theory. In this subsection we argue how the results 
discussed above apply to some additional cases listed in table [2j 

The potentials Vq and Vj i have been introduced in The former was used to mimic 
the QCD equation of state using holography while the latter have a second-order phase 
transition at T = T c . Fixing b = 0.606 corresponds to = 0.15 in the infrared [12] . 
The variety of potentials leads to a range of conformal weights 3 < A < 3.93 (table [2]). 
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potential 

b 

c 2 

c 4 

c 6 

A 

^6 

0.606 

2.06 

0 

0 

3.93 

v 7 

0.606 

1.8 

0 

0 

3.79 

bs 

0.606 

1.145 

0 

0 

3.37 

bg 

0.606 

0.808 

0 

0 

3.10 

Vio 

0.606 

0.703 

0 

0 

3.00 

V n 

1/V2 

1.942 

0 

0 

3.37 


Table 2: The second group of potentials considered in the present paper. 

r 2 Imw 




Figure 5: The speed of sound c 2 (left) and the imaginary parts of the lowest quasinormal 
modes at k = 0 (right) for the potentials from table 2. The one with cusp in the c 2 
(left) corresponds to the most decreased in damping (right). 


For all these cases the qualitative conclusions discussed earlier in this section still hold. 
In figure [5] we plot the speed of sound (left) and the imaginary parts of the lowest 
damped QNM’s (right) for the potentials Vq — V\\. Note that in all cases there is a 
critical temperature, corresponds to the lowest value of c 2 , which might be considered 
as the cross over/phase transition point (related to the potentials). Damping of the 
lowest QNM’s decreases close to the T = T c by a factor of 2 — 7 (depending on the 
potential) relative to the conformal theory at high temperatures, figure [5] (right). 

Figure [6] shows the momentum dependence of the imaginary part of QNM’s (left) and 
the ratio of the real parts to the imaginary parts (right). The former indicates that 
ultralocality still holds for this more diverse class of potentials. Interestingly, comparing 
the plots in figure [5] (left) and [6] (left) suggest that the phenomenological relation we 
found in previous section^] namely 


z 2 ~ 1.114 


Im u"(k = 0) 
2v tT 


( 36 ) 


is valid for this class of potentials too (again untill T T c ). 

5 The coefficient 1.114 is the same as ”a” introduced in figure |4j 
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Imoi" (ft = 0) 


2nT 



| Imw | 



Figure 6: The curvature at k = 0 for the damping frequency for potentials from table 2 
(left) and the ratio of the real to the imaginary part (right) for the same potentials. 


The ratio of real parts to the imaginary parts of the lowest QNM’s in figure [6] (right) 
shows a decreasing at T ~ T c i.e. the lower c 2 at T = T c the bigger decreasing in 
Rew/Imw. 

Surprisingly, even the potential Vn with the second-order phase transition exhibits 
qualitatively the same behaviour as the other potentials discussed in the present paper. 


6 Conclusions 

In this paper we carried out a study of the lowest quasi-normal modes in a class of 
nonconformal holographic models which exhibit an equation of state very similar to 
the one obtained using lattice QCD. This class of models was introduced in ra and 
incorporates 5-dimensional gravity coupled to a scalar field with a given self-interaction 
potential which parametrizes the model. 

The frequencies of the lowest quasinormal modes provide a scale for the importance 
of nonhydrodynamic degrees of freedom, thus their determination is of a definite phe¬ 
nomenological interest. Our main observations are the following. 

Firstly we found that, within the class of considered models, the imaginary part of 
the QNM frequency is strongly correlated with the speed of sound characteristic of the 
equation of state, once we factor out the trivial conformal temperature dependence. In 
particular, it decreases by a factor of around two at the point of the QCD crossover 
transition. This means that nonequilibrium effects will become more pronounced closer 
to the QCD phase transition/crossover. This seems to be a robust characteristic of this 
class of models and persists also for models with other equations of state considered 
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for completeness in section 5.4 We provided a phenomenological formula (33) linking 
the damping with the speed of sound. It is important to emphasize, that although the 


numerical values of the coefficients in (33) are specific only to the models mimicking 
the QCD equation of state, similar fits, with coefficients of the same sign and similar 


order of magnitude, work also for other models considered in section 5.4 


Secondly, we found that the quasinormal frequencies practically coincide for a whole 
class of models (potentials) which lead to the same equation of state (or more precisely 
to the same speed of sound c 2 [T) as a function of temperature). In particular, they 
seem to be quite independent of the particular UV properties of the concrete potential 
such as the anomalous dimension of the operator deforming the theory. 


Thirdly, we found that the property of ultralocality found in CD, namely the very 
mild dependence of the quasinormal modes on the momentum k around k — 0 persists 
away from conformality. Even more so, it becomes more pronounced. We also noticed 
an intriguing feature that the curvature of the imaginary part of the QNM frequencies 
around k = 0 follows surprisingly well the speed of sound squared c 2 . 


We believe that the above observations should be of phenomenological interest, espe¬ 
cially as they indicate a more pronounced role of nonhydrodynamic degrees of freedom 
close to the QCD phase transition/crossover. It would be very interesting to gain some 
analytical understanding of these properties as well as to investigate directly nonlinear 
dynamical evolution in such models. 
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